










* 'vc V .I<i'.'i''“4'%‘’ 




W&P? 


't, '*/-'■ '/if. 


•i. ;"-:t;i?L-..l-A 


RESPONSE OF A NONLINEAR AIRCRAFT TO RUNWAY INDUCED 
EXCITATIONS DURING VARIABLE VELOCITY RUNS 


A Thesis Submitted 

In Partial Fulfilment of the Pequirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

KAPADIA KAIKOBAD E. 


to the 

DEPARTMENT OF AEROSPACE ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

MAY, 1989 



^ 3 OCT W89 



Acc. No- 






KAP- 




This is to certify that the work 'Response of a Nonlinear 
Aircraft To Runway Induced Excitations During Variable Velocity 
Runs’ has been carried out under my supervision and has not been 
submitted elsewhere for a degree. 


(DR. DAYANAND YADAV) 
Assistant Professor 
Department of Aerospace Engineering 
Indian Institute of Technology, Kanpur. 




Acknowledgement 


I wish to express my deep sense of gratitude and thankfulness 
to Dr. Dayanand Yadav for his sustained efforts in the 
completion of this work in the present form. During the entire 
tenure of this work he has been a constant source of inspiration, 
guidance and encouragement to me. 

I am grateful to my hall residents who made my stay at 
I. I. T. Kanpur a very memorable and pleasant one and above all to my 
friends Ramamoorthy, Gautam, Sanjay and Bohra in particular, who 
have helped me at various stages of my work. 

I am thankful to Mr. A. K. Ganguly for his painstaking efforts 
in preparing the drawings and to Mr. G.R. Boshing for his excellent 
typing of the manuscript . 

Thanks are due to everybody in my department who provided a 
very warm and friendly enviornment to carry out this work. 


(KAPADIA KAIKOBAD E. ) 


(ii) 



Abstract 


Response of an aircraft in variable velocity runs during 
landing, take off and taxiing on runways, with non~aero mean and 
random unevenness, is studied. The aircraft is modelled as a rigid 
body and the runway is considered as a rigid pavement. The study 
incorporates the mathematical modelling of the aircraft including 
the landing gear and the runway. 

The modelling of the landing gear includes non-linearity in 
the spring and the damper. It incorporates the polytropic 
compression of the air, compressibility of the hydraulic fluid and 
the internal frictional forces in the shock struts. The nonlinear 
tyre spring and the aerodynamic lift are also accounted for. The 
effect of inclination of the shock struts and the linkage dynamics 
is included in the analysis. For the present study, the model is 

limited to a rigid heave model with single point input from the 
track. 

The track unevenness is prescribed as a random process in 
spatial coordinates over a mean level. The random unevenness is 
described by its power spectral density (p.s.d.). The mean level 
may include slopes, bumps, dips. ramps, bomb crater repair patches, 
step differences, etc., singly or in combination to simulate 
different possible faults in the runway. 


Monte-Carlo simulation approach is used for the solution of 
the problem. The system dynamic equations have the initial 
conditions prescribed by the physical parameters . The track input 
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forms the forcing function. The input consists of random unevenness 
superimposed over deterministic mean. 

The equations of motion along with the forcing function are 
integrated numerically to find the response of the system. 

Response ensemble is generated by repeated solutions with 
different track input samples. The second order statistics of 
response is obtained by averaging across the sample sets in the 
ensemble. 

Numerical data is presented for landing and take-off runs with 
parametric variation. The parameters studied are the roughness 

level of the track for landing and take-off and the sink velocity 
in landing. 
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Nomenclature 


length of the drag linkage AE. 

length of the part of the main linkage carrying the main 
shock strut. 

length of the entire main linkage. 

distance between the two points of the linkage hinged to 
the aircraft frame. 

perpendicular projection of the main linkage. 

pneumatic area. 

hydraulic area. 

orifice area. 

roughness constant. 

coefficient of discharge. 

pneumatic force. 

journal friction force. 

seal friction force. 

hydraulic force. 

horizontal force on the wheel axle. 

horizontal ground reaction. 

axial force in the drag shock strut. 

axial force in the main shock strut. 



preloading of the drag strut. 



F preloading of the main strut. 

s 

20 

F vertical force on the wheel axle. 

V 

a 



vertical ground reaction. 


L 

I 
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n 



lift force acting on the sprung mass. 

axial distance between the shock strut lower bearing and 
the lower end pin joint of the piston tube, 
axial distance between the upper and lower bearings on 
the extended shock strut. 

polytropic exponent for the air compression process, 
initial air pressure 



stroke of the drag strut, 
stroke of the main strut, 
initial air volume. 



initial oil volume. 


W 


W_ 


Z(t) 


weight of the sprung mass, 
weight of the unsprung mass. 

horizontal relative displacement of the sprung and 
unsprung mass, 

Gaussian Random state vector. 


sprung and 


vertical displacement of the sprung mass, 
vertical displacement of the unsprung mass, 
relative vertical displacement of the 
unsprung mass, 
oil bulk modules. 

angular position of the main linkage at any instance of 
time. 

initial angular position of the main linkage, 
coefficient of friction of the track, 
coefficient of friction of upper bearing. 

coefficient of friction of lower bearing. 

coefficient of seal friction, 
mean of a random process x. 

correlation constant. ' 

mass density of the hydraulic fluid, 
standard deviation of a random process x. 
power spectral density. 

angular position of drag linkage at any instance of 
time. 

initial angular position of the drag linkage, 
spatial frequency. 

Gaussian sequence, 

differentiation with respect to time. 
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Chapter I 


Introduction 

1.1 Problem. 

All conventional aircrafts spend a significantly large part 
of their useful life during taxiing, take-off and landing on 
runways. During these traversals on runways, the aircraft follows 
the ground profile. The ground unevenness, therefore acts as a 
base excitation to the aircraft during these motions. Besides the 
nature of the ground unevenness, the induced excitation also 
depends on the type of contact and the traversal velocity of the 
aircraft. 

The vibration of the aircraft during its operation on the 
runway induces dynamic loads on the structure. This has become 
more and more important due to the increased need of reducing the 
time spent by an aircraft during landing, take-off and taxiing. As 
a consequence of this the aircraft traversal velocities have 
increased. These vibrational loads cannot be ignored, since their 
cumulative effects represent a considerable proportion of the 
total fatigue damage. The fatigue life of an undercarriage at its 
points of attachment with the airframe, and the associated 
structural components make up for an important design case. 

Apart from the metal fatigue there are the equally important 
dynamic peak stresses that arise in a fully laden aircraft during 
its manoeuvres on rough runways. Thus it is important for 
designers to be aware of the vibrations induced due to runway 
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unevenness, so that an attempt can be made to reduce them. 

The effects of these vibrations on passengers and the crew 
is very important. It is known [1] that cockpit vibration levels 
are roughly twice than those experienced at the aircraft mass 
centre, and thus pilots are more influenced than the passengers. 
In some cases high cockpit vibration levels in the vertical plane, 
can degrade the pilot’s ability to read instrument panels 
accurately, and apply the correct amount of control required. It 
is also very important from the point of view of passenger 
comfort. 

The dynamic displacements react through the oleo -pneumatic 
units and tyres with the runways which induce vibrations, thereby 
creating additional dynamic loading of runway surfaces. This 
results in their general deterioration, which would in turn again 
lead to an increased level of vibration causing a viscious 
derogj^atory cycle. 

The rough ground can be treated either as a homogeneous 
(spatially stationary) or as a non-homogeneous random process 
depending on its statistical description. For the runway 
unevenness described by a stationary random process, the response 
of the aircraft can be predicted by classical stationary random 
vibration analysis, if the aircraft is moving at a constant 
velocity, as is the case commonly during taxiing. Contrary to 
this, during landing and take-off, when the aircraft accelerates 
or decelerates, the excitation becomes non-stationary even though 
the runway unevenness is stationary. Non-stationary excitation may 
also result with non-stationary ground profile. 
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For an accurate modelling the problem should take into 
account the following features : - 

(a) nonhomogenity of the ground profile. 

(b) nonlinearity of the aircraft dynamic characteristics. 

(c) the aircraft traversal at both constant and variable 
velocities. 

1.2 Review o/ the Past Worh. 

This section briefly outlines work done by other researchers 
in the mathematical modelling of the landing gear and different 
approaches used for obtaining the response of the aircraft to 
runway induced excitations . 

1.2.1 Mathematical Modelling 

As early as 1948 Stowel, et. al . £2] investigated the effect 
of flexibility of the wing on the landing impact forces. The 
aircraft was modelled as a simplified structure consisting of a 
uniform bar for the wing, a concentrated mass for the fuselage and 
an undamped linear spring for the landing gear. It was found that 
for moderately flexible landing gears neglecting structural 
flexibility led to small conservative errors. 

Later Cook and Milwitzky [3] studied the problem by 
considering several free-free modes of vibration of the aircraft 
and the non-linear characteristics of the landing gear. They 
indicated that the errors involved in neglecting the flexibility 
are not conservative. 

Raghavan [43 also concluded that the maximum response after 
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landing impact is 5 to 29% less for a flexible system. 

In view of the above mentioned investigations, it appears 
that the landing gear performance is relatively unaffected by the 
elastic deformations of the airplane structure. In addition, the 
proximity of the landing gear to the nodal points makes the 
conclusion all the more acceptable. 

Milwitzky and Cook [5], probably for the first time, 
investigated the behaviour of the conventional type of 
oleo-pneumatic landing gear during landing impact. The basic 
analysis is presented in a general form and motions of the landing 
gear prior to and subsequent to the shock strut deflection were 
taken into account. The first phase of the impact was considered 
as a single degree of freedom system in order to find the initial 
conditions of motion at the instant of shock strut deflection, 
after which the landing gear was considered as a two degree of 
freedom system. The effects of the non-linearities of the landing 
gear parameters and that of the drag forces on the wheel were 
included. It was concluded that the orifice discharge coefficient 
has a marked effect on the behaviour of the landing gear, whereas 
the air compression was found to be relatively unimportant. Linear 
tyre characteristics gave good results during normal impact. 

Refinements were added by Wahi [6,73 by including orifice 
function, hydraulic fluid compressibility and polytropic exponent 
of air during the compression stroke. 

^ ^ ^ ^ ^ Reddy, et. al. [83 carried out analysis on a 
semi-articulated gear with geometry included as a part of the 
analysis. They also included oil compressibility effects. To 
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identify the active parameters, they carried out a parametric 
study of various parameters like coefficient of friction of main 
orifice C^, polytropic exponent of the air compression process n, 
coefficient of friction at lower and upper bearings /j . & /j . 

etc. The ground reaction factor Fj^ defined as the ratio of the 
ground reaction force to the wieght of the aircraft to be 

e 

supported by the landing gear W was taken as the objective 
function. The travel of the wheel axle was restricted to a value 
of 255 mm during landing. 

It was shown that there is a value of C^, for which the gound 
reaction factor is minimum. The ground reaction factor is not 
significantly affected by a variation in the polytropic exponent. 
The internal friction force in the shock strut remains nearly 
constant with changes in coefficient of friction. 

1.2.2 Response Analysis 

Traditionally undercarriages have been designed to absorb 
landing impact forces rather than vibratory forces caused by 
runway unevenness. For a relatively small aircraft, runway 
excited vibrations were considered unimportant. Thus despite the 
growth of the aircraft size little attempt has been made to design 
undercarriages from the vibration point of view, with the result 
that vibration levels have increased with aircraft size. Not much 
analytical work exists in the literature, may be because of the 
fact that the non-linear random vibration theory had not reached a 
stage SO as to provide a tractable approach to consider inherent 
non-linearities of the aircraft. 
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It was only in the early sixties that Silsby [10] studied 
analytically the effects of aircraft and landing gear parameters 
on the airplane response to runway roughness. He considered a 
highly simplified linear model for the landing gear and a rigid 
airplane structure. The results of the study indicated that, 
within the limitations of the assumptions made, the parameter 
variations had little effect on the motion of the centre of 
gravity of the aircraft, whereas the motion of the pilots’ 
compartment was strongly influenced by the speed. 

Tung, et. al. [11] have computed the dynamic response of 
supersonic transport aircrafts to runway roughness during taxiing 
and take-off through a deterministic approach, in which the 
response is computed by numerically integrating the equations of 
motion of free-free aircraft. Nonlinearity of spring and damping 
forces have been taken into account. A statistical approach is 
also highlighted in which the runway unevenness is represented as 
a stationary random process. The major obstacle to applying a 
statistical analysis was that of the nonlinearities present in the 
oleo- pneumatic undercarriage, in the form of air spring stiffness, 
velocity squared damping and Coloumb friction. It was not 
possible to directly apply the theory to nonlinear systems and so 
two important approximate techniques are shown to apply. The 
methods are : 

(a) The Perturbation Technique developed by S.H. Crandall. 

(b) The equivalent linearization technique developed by 
T.K.Caughey. 

Both mehtods assume that the input is a stationary random 
process and the nonlinearities are small. The perturbation theory 
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can only be applied ■bo nonlineari'ties which can be represen'ted as 
polynomials e.g. , the air spring. But since the V|V| damping and 
the Coloumb friction cannot be represented by polynomials, 
perturbation theory cannot be thus employed. So it is suggested 
that first apply the linearization technique to obtain the 
equivalent linear damping coefficient for the nonlinear damping 
forces for the aircraft, considered as a single degree of freedom 
system, and then apply the perturbation method to a multi-degree 
of freedom system with an equivalent viscous damping coefficient, 
containing a polynomial representation for the nonlinear air 
spring force. The means for implementing such a scheme, however, 
was not given. But the paper applies the perturbation method 
to a two degree freedom rigid heave model. The nonlinear air 
spring was represented by a polynomial but only viscous damping 
was assumed. The paper concluded that although the equivalent 
linearization method was suitable for a single degree of freedom 
system with nonlinear damping, it appeared unlikely that it could 
be applied to a multi-degree of freedom system. 

Kirk and Perry [12] have carried out the analysis of taxiing 
induced vibrations of a subsonic aircraft by the power spectral 
density (p.s.d.) method. They considered one elastic mode coupled 
with the rigid mode of the aircraft. They converted the nonlinear 
damping into an equivalent linear damping. In using this, 
approach, one of the basic assumptions which is open to question, 
eoncerns the stationarity of runway profiles. Evidence regarding 
stationarity of runway profiles is given in reference £12], for 
relatively smooth surfaces (i.e., good quality runways in 
commercial use) having rms displacements less than about 0.018 m. 



in which case the probability density function can be considered 
as Gaussian. For rougher surfaces, however, Gaussian distribution 
could not be assured without further investigations. 

When the aircraft has a variable velocity, as is the case 
during landing and take-off runs the input excitation becomes 
nonstationary. Virchis, et. al. [13] and Sobczyk, et. al. [14] 
have used the time domain approach to determine the response 
statistics for a single degree of freedom model in constant 
velocity and constant acceleration runs. 

Yadav and Nigam [15] use frequency domain formulation for 
response evaluation . The equations of motion are transformed to 
space coordinates, resulting in the governing equations having 
space dependent coefficients. The input process is spatially 
stationary. Due to this the solution admits an evolutionary 
spectral form of representation and for a particular velocity 
profile (viz., S{t) = a^^ + closed form solutions for 
response statistics are presented after some lengthy integration. 
Numerical mehtods could again be employed to handle more general 
situations. 

Hammond, et. al. [16] computed the non-stationary response of 
a vehicle modelled by linear dynamics travelling over a 
homogeneous gound profile at variable velocities. The method uses 
linear state space technique and incorporates one important 
assumption, that the gound surface undulation may be described as 
the output of a white noise excited shaping filter (in the spatial 
domain). A set of first order equations are integrated 
numerically to obtain the zero lag response covariances 



9 


(correlations) , 

Later in another paper Harrison and Hammond [17] obtained the 
time varying mean and covariance of nonlinear systems excited by 
non- stationary random process, the non-stationarity arising out of 
the variable velocity of the vehicle. They used statistical 
linearization to overcome the nonlinear problem. 

1 . 3 Fr&sent Work 

For the present work, the response of an aircraft in variable 
velocity runs, as in taxiing, take-off and landing, is studied. 

For a simplistic analysis, the weight of the aircraft to be 
carried by the landing gear is considered as a concentrated mass 
and the weight of the wheel unit is also assumed to be 
concentrated at the wheel axle. The input from the track 
unevenness is considered to be a single point input as only one 
landing gear is studied. These assumptions limit the aircraft 
model to a rigid heave model with a single point input from the 
track . 

The landing gear is modelled as a telescopic landing gear 
with two shock struts, one the main strut and the other the drag 
strut. To make the model more general both, the main and the drag 
struts are inclined with the vertical and the geometry of the 
linkage is also included. The modelling includes nonlinearities in 
the spring and the damper. It also incorporates the polytropic 
compression of air, compressibility of hydraulic fluid and the 
internal friction forces in the wheel struts. Nonlinearity in the 
tyre spring and the aerodynamic lift is also accounted for. 



so that 


The runway track Is assumed to be a rigid pavement, 
in the absence of any flexibilities it does not interact 
dynamically with the response analysis. 

The ground unevenness is prescribed as a random process In 
spatial coordinates superimposed over a mean level. The mean level 
includes slopes, bumps, dips, bomb crater repair patches, ramps 
etc. , representative of the faults in the runway. 

The random track profile is numerically simulated to fit the 
track description in the following steps 

(a) Initially, an uniformally distributed, non- repeatable 
sequence of random numbers between zero and one is 

(b) A Gaussian white sequence, having zero mean and unit 
variance is generated by using the above numbers. 

(c) This Gaussian sequence is operated upon by two 

differential operators in sequence to produce the 
required Gaussian random process. 

After the simulatlou of the Gaussian process from the 
prescribed power spectral density (p.s.d.), a numerical scheme is 
used to obtain the p.s.d. of the simulated process to verify its 
closeness to the target. Monte-Carlo simulation approach is used 
for the solution of the problem. The system dynamic equations have 
the initial conditions prescribed by the physical parameters. The 
track input forms the forcing function. The equations of motion 
along with the forcing function are integrated numerically to find 
the response of the system. The response samples are generated to 

form an ensemble and then analysed to determine its statistical 
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parameters . 

During the course of this work, response samples have been 
generated for analysis during landing and take-off, from a runway 
having a constant mean slope. For the study of the landing case, 
variation in the roughness of the runway track and variation in 
the sink velocity has been considered. For the take-off case 
variation in the roughess of the runway profile is studied. > 

1.4 Strxictxjora.1 Layout of the Thesis 

The second chapter embodies the system modelling, modelling 
of the runway track and the system equations of motion. In the 
third chapter, the approach used for the numerical solution is 
discussed. The fourth chapter constitutes the result and 
discussion part of the thesis and the fifth chapter concludes the 
present work with some suggestions for futuristic work. 



Chapter II 


Theory 


2.1 Introduction 

VHienever an aircraft lands, takes-off or taxiies on the 
runway, its dynamics response is mainly affected by the various 
landing gear parameters, because the ground input acting as a 
forcing function is interacted by the landing gear system. 


In the 

present work 

only 

the vertical 

heave motion 

is 

considered; 

The combination of 

airplane 

and 

landing gear 

is 

idealised to 

constitute a 

two 

degree of 

freedom system. 

The 


landing gear used is of the articulated type with oleo-pneumatic 
shock struts. The equations of motion are developed to include 
the non-linearities of the landing gear parameters, the wing lift 
and the inclinations of the landing gear linkages carrying the 
oleo-pneumatic struts. In the analysis, the motion of the landing 
gear prior to and subsequent to the beginning of stock strut 
deflection is treated. Monte-Carlo simulation approach is used for 
generating the ground input. The solution for the equations of 
motion along with the input, is obtained by employing a numerical 
scheme of integration. 

2.2 System. Modelling 

2.2.1 Representation of the aircraft 

For an initial analysis, the airplane is considered as a 
rigid lumped mass with only heave degree of freedom. As the 
airplane is constrained to move only in the vertical direction, a 
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single landing gear with its share of the aircraft weight W1 can 
be considered for the purpose of the analysis. Thus the landing 
gear together with the airplane mass constitute a system shown in 
Figure 2.1. Assuming the weight of the lower portion of the 
landing gear and wheel, W2, to be concentrated at the wheel axle, 
the airplane reduces to a two degree of freedom system. 

The present analysis is subjected to the following 
assumptions : 

i) Airplane is a rigid mass. 

ii) Mass of the landing gear and the wheel with its 

attachment or the unsprung mass, is concentrated at the 
wheel axle . 

iii) Shock struts are infinitely rigid in bending. 

iv) Only vertical motion of the aircraft is considered. 

2.2.2 Mechanics of the Gear 

Figure 2.1 shows the system model used for analysis in the 
present work. It consists of two linkages AE and BD carrying the 
drag and the main oleo-pneumatic struts. The drag linkage is 
hinged at both ends, one end to the aircraft frame. A; and the 
other to the hinge point E on the main linkage. The main linkage 
has its upper end B hinged to the aircraft frame and carries the 
wheel set at its lower end D. The arm CE is rigidly connected 
perpendicular to the main link. A^ to A^ are the geometrical 
parameters describing the gear. 

The vertical travel of the aircraft is and the vertical 
deflection of the tyre ie Z,, both measured positive 




downwards from the unloaded fully extended condition of the gear 


2.2.3 Kinematics of the Gear 

The kinematics of the telescopic gear is shown in Figure 2.2. 
At any time instant t, the linkages AE and BD move from their 
original position to the deflected position AE’ and BD’ . Z and 
x^ represent the differential travel of the wheel axle with 
respect to the main mass in the vertical and the horizontal 
directions respectively. 

The kinematic relations for x^ , , main strut closer, 

and their derivatives are established as presented below. 

From geometry - 



tt 

. cos & 

o 

A 

11 

. sin & 

o 

A 

<5 

fl 

. cos & 

o 

As 

11 

. cos <t) 

o 


ft 

. sin 

o 


( 2 . 1 ) 


As Z^ is the differential movement at the wheel axle point 


Z . = Z, - Z 
a 2 1 

z, = z. - z 


z^ = z_ - z 


( 2 . 2 ) 


Also, let the new position of the wheel axle in the vertical 
plane be given by 
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Then the shock strut stroke and the strut telescoping 

velocity are given as ; 


i/z 

CX^ + + u' j <2. 4) 


( X , +A_ ) X . + uu 

• — d 15 a 

and the acceleration of the strut S can be written as 

z 


1/Z 


(2.5) 


Cx,+A^)x. + uu 

a 5 a 


s. = 




cx^tA^r t 


“1 


3 /Z 


Cx^+A^)x^ + (x^)* + uu + (u)* 


Cx^+A^)‘ -H 




1/2 


( 2 . 6 ) 


From the gear geometry the angle of inclination of the main 
struct (&) is expressed as 

® ^ ] 

Also from the figure 2.2 it can be seen that 

[ -| 1 / z 
CA,-Sp^ - ] 

sin O = 

(^3-5^) 

Differentiating equation (2.7) and substituting the value of 


(2.7) 


( 2 . 8 ) 


sin 0 from equation (2.8) gives 
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6 = 


- uS^- u (A 3 -S^) 




1/2 


(2.9) 


Differentiating equation (2.9) once more 7 ields 


.. - uS - u (A -S ) 

e - = = — = — 


SJuS^+a(A3-S^ )] 


(A3-S^)[ (A3-S^)"-u"]‘"" 


(A3-S^)^[(A3-S^)"-u=']""" 


[ (Ag-S^ )S^+uu3 [ uS^+u(A3-S^)] 

’ (A3-S^ )[ (A3-S^)"-u"]""^ 


( 2 . 10 ) 


The horizontal displacement of the wheel axle x. is 

Cl 


X, = O’D’ - OD 

d 

= (A3-S^)sin© 
= (A3-S^)sin e 


A^sin 

3 o 


- A 


( 2 . 11 ) 


On differentiating it once, the velocity x, is 

Cl 


X. = (A^-S, )cos Q .O - S, sin 0 

O 3 Z Z 

= u© - S,sin e 
2 


( 2 . 12 ) 


On differentiating euation (2.12), the corresponding 
acceleration x. is 

a 


X . = u© + u© - S^cos ©.© - S,sin © 
d 2 2 


(2.13) 
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The kinematic relations for the stroke of the drag strut S 
and its derivative is established as follows. 


Looking at figure 2.2 it can be seen that 


BD’ = A = A -S 

3 1 3 2 

BG’ = A = A -S 

2 1 Z 2 

0’ ’ ’C’ = BC^ . sin e 


= (A^-S^ ) . sin e 


BO’ ’ ’ 

AP’ = to - (A^^ .sin 0 ) 

E’P* = A^ - 0”’G> - (A^^.cos e) 


= to = ^^2~^Z ^ 


A 


= [ AP’* + E’P’*] 


j - y-y 3 ^ *1 1 ^ ^ 


(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 


••• S, = A^-A t 


Also tan <p = E’P’/AP^ 


( 2 . 21 ) 

( 2 . 22 ) 


Differentiating equations (2.16), (2.17), (2.18) and (2.19) 
once each we have the following derivatives 


O’ ”C’ 


(A^ -S^ ) . cos 0 .<9 - (S^sin©) 


= -(A,-S, )sin & .& - (S cos 9) 

Z Z 2 

AP’ = to - (A^^.cos 0.(9) 

lO 

E’P’ = -0”’’C’ + (A^^. sin 0.0) 


(2.23) 

(2.24) 

(2.25) 

(2.26) 


Differentiation of equations (2.21) and (2.20) subsequently 
yields 


S = - A = 

1 It 


- C(AP’ . AP’) + (E’P’ . E’P’)) /A 


It 


(2.27) 
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These expressions are used for calculating the shock strut 
forces required to solve the equations of motion. 

2.2.4 Forces in the Shock Strut 

It can be observed that the system under the present study 

consists of two oleo -pneumatic shock struts. One strut is located 

on link BC and is the main strut, the other one located on the 

link AE is the drag strut. During the compression stroke the 

hydraulic fluid is forced out of the hydraulic chamber through an 

orifice to a lower chamber, providing a nonlinear damping which is 

.2 

proportional to the square of the strut telescoping velocity S. 
When the fluid enters the lower chamber, it pushes a separator 
downwards in the pneumatic chamaber, thereby providing a nonlinear 
spring due to the polytropic air compression in the pneumatic 
chamber. 

The shock strut axial force is a combination of the pneumatic 

force F , the hydraulic force F. , the journal friction force 
a n 

F- and the seal friction force F- [8] and is expressed as 
j s 

F = F + F, + F„ + F„ (2.28) 

s a xi X • I 

J s 

The pneumatic force is 

a 
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n 

(2.29) 


where 

v„ = initial air volume 

a 

o 

= initial oil volume 

o 

o 

P„ = initial air pressure 

^o 

P = instantaneous air pressure 

Ql 

ft = oil bulk modulas 

n = polytropic exponent for the air compression 

process . 

while computing the pneumatic force F , the compressibility 

s. 

of oil has also been considered. [7]. The solution of the 
transcendental equation (2.29) is obtained by the method of 
successive approximation. 

The hydraulic force F is given as [53 

3 . 



|S| 2(C^A^)* 

where 

p = mass density of the hydraulic fluid 
= hydraulic area 
= coefficient of discharge 


(2.30) 


orifice area 
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The Journal fric'tion force in hhe s'tru't xe *the sum of 'the 

3 

journal friction forces contributed by each of the two bearings . 
It is obtained as [5] 


where 




+ P . ) 
"^2 





(2.31) 


F 


n. 





normal load on the shock strut 

coefficient of friction of upper and lower bearings on 
the shock strut 


= axial distance between the upper and lower bearings on 
the extended shock strut 

= axial distance between the shock strut lower bearing 
and the lower end pin joint of the piston tube. 


A consideration of the system shown in figure 2.1 will show 
that the drag strut on the link AE is always axially loaded and so 
this strut, the journal friction force does not come into 
picture. 


Also, it can be observed that the normal force on the main 
shock strut is given as 

= F^ cos© sin© (2.32) 


Substituting equation (2.32) in (2.31), the journal friction 


force is 
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(2.33) 


where 



S 

|S| 


I cos 6 . sin 6 1 


3 ^ 3 


) 

2 


^ +s 

1 


+ ^ 



The seal friction force is assumed to be proportional 

s 

to the pneumatic force [8] as 



(2.34) 


where 

fj = coefficient of seal fx'iction 
s 

By substituting equation (2.33) in equation (2.28) and 
rearranging, the shock strut axial force is obtained as 


= (F„ + ) / 


a 


n 




(2.35) 


2.2.5 Forces on Tyre 

The force transmitted from the ground to the tyre, F can 

g 

be expressed as a polynomial in tyre deflection with the 
coefficients obtained to fit the experimental curve provided by 
the tyre manufacturers. 
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s 

+ a {Z +h)'‘ + a^{Z +h)° + a^(Z,+h)® 

+ a (Z +h)’’ + a„(Z +h)® + a„{Z,+h)'’ 

+ (2.36) 

where, 

h = profile of -the ground surface 

= displacement of the unsprung mass 
a^ , , . . . , = constants 

The instantaneous drag force between the tyre and the ground 
F, is assumed to be equal to the vertical ground reaction F 
multiplied by a suitable coefficient of ground friction p 

®'h = • ®'v ^2.37) 

g ^ g 

Thus the vertical and horizontal forces on the tyre at the 
place of ground contact are determined from equations (2.36) and 
(2.37). 


2.2.6 Equations of Motion 

The dynamic system considered in the analysis is shown in 
figure 2.3. The portion of the aircraft weight acting on the gear 
or the sprung mass is W1 . The various forces acting on it in the 
vertical direction are shown in this figure. The mass of the 
shock strut pistons and linkages is transferred to the unsprung 
mass, which also consists of the mass of the wheel unit assembly. 



s various forces acting on the unsprung mass are also shovm in 
s figure 2.3. 

The shock struts on both the main and drag linkages have 

sloads F and F respectively due to the extended or initial 
°20 ®10 

c pressures pa^ and pa respectively. Until the preload of the 
z 1 

rut is exceeded the strut acts as a rigid member. 

Initially, until the preloads of the two struts are exceeded, 
Ly the tyre deflects, and the system acts as a single degree of 
sedom system with = Z. The equilibrium equation of the 

stem is 

•|z + L- W + F =0 (2.38) 

g 

After the preload in either the main or the drag shock strut 
exceeded, the system starts behaving as a two degree of freedom 
item with either the main or the drag or both the struts 
Lescoping. The vertical force equilibrium of the unsprung mass 

^ + F - F - = 0 (2.39) 

g a 

I the vertical force equilibrium of the sprung mass is 

^ \ + L - W^+ (F^ .cos <p) + (F^ .cos e) = 0 (2.40) 

The horizontal equilibrium of the unsprung mass gives 
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0 


( 2 . 41 ) 


where, 

= horizontal force on the wheel axle, 
a 

2.2.7 Initial Conditions 

To establish the point at which the transition from one 
degree of freedom system to the two degree of freedom system takes 
place, it is necessary to solve equation (2.38) for the single 
degree of freedom system. 


Equation (2.38) can be rewritten as 


Z = (W-L-F^ ) / (W/g) 
g 


( 2 . 42 ) 


The vertical force equilibrium of the linkage system can be 
written as 


F = F cos <^ + F cos & 
vs s 

a 1 2 


(2.43) 


Substituting equations (2.42) and (2.43) in (2.40) yields 


F = F 
a g 


( 2 . 44 ) 


The horizontal force equilibrium of the linkage gives 


sin © + F, ” F„ sin 0 = 0 
s n s 

2 a 1 


( 2 . 45 ) 
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Equation (2.41), after substituting equation (2,37) and 
neglecting the inertia force acting on the unsprung mass, yields 

F, = pi F (2.46) 

Substituting equation (2.46) in (2,45) gives a relationship 

between F and F^ as 
s s _ 

1 2 

F = (F sin ?!» - F„ )/sin © (2.47) 

S S g V 

2 1 g 

Substituting equation (2.47) in (2.43) we have 
CF^ +(pig F^ )cos ©/sin ©) 

F = 5 i (2.48) 

The calculated values of F and F, in the one degree case 

s s 

2 1 

are compared with their corresponding preloads at each step. The 
instance at which any one of these exceeds the preloading 
transition takes place from the one degree system case to the two 
degree case. 

2.3 Modelling of the Runujay Profile 

As has been stated before, the input to the landing gear 
system arises due to the motion of the aircraft on the runway, 
which inherently has roughness to a greater or a lesser extent. 
It is this undulating ground profile which acts as a forcing 
function for the undercarriage system. During the present work, 
the ground profile is numerically simulated as a Nonstationary 
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Gaussian random process. The profile h(s) in space coordinates is 
expressed as the sum of a deterministic and a random part. 

The deterministic part may arise due to a variation in the 
mean of the track level because of faults such as aclivity, 
declivity, bumps, dips, step-ups, step-downs, ramps and repair 
patches resulting due to improper leveling, differential 
settlement of slabs, repair at joints and bomb craters. Some of 
these faults are shown in figure 2.4 and can be used singly or in 
combination to model the various faults in the runway. 

The random profile generated represents the ground 
undulations. The runway unevenness is considered to be a Gaussian 
Random process with a power spectral density as 


§ia) 


Ci^ + 


( 2 . 49 ) 


where , 

^ = spatial frequency 
C = roughness constant 
V = correlation constant 

The algebraic sum of the determinstic profile and the random 
profile gives the non stationary runway profile. 

2.4 Response Statistics 

Repeated solution of the system equations can be resorted to 

The second order 


generate an ensemble of the system response. 


30 


I 



cuT 


r-rr-rr-rT-rT-T 


CONSTANT MEAN SLOPE 



DIP IN LEVEL 



A BUMP 


k- ^ ' .1 

\ 

r n 

7"7y/v 

oC 

V ; H 

V v r; / / / 


X 




n 

7T~T~7~7~7~r~T7~ 

H 

Ti / i / i J ? / / / / /-^ 

X 



7 ^ ^ / /' / / / / / / r 

/ 


r 

H 

[ 


' ^ /////// / 


REPAIR PATCH 


STEP UP 


STEP DOWN 


FIG. 2.4 FAULTS IN THE RUNWAY TRACK 


31 


s'ba'tis'tics of "the Tesponse can "bh-en be obtained, by averaging 
across the ensemble as [19] 

1 N 

1 E x^(t) 

1=1 

2 1 N 

°'x = iFT ^ CX^(-t) - 

1 = 1 

where, .^i^and denote the mean and the standard deviation 
of a random process X with N as the number of samples in the 
ensemble. 


(2.50) 

(2.51) 



Chapter III 


Numerical Solution Approach 

3.1 IntrodvLCtion 

To ob-tain the response of the aircraft during its runs on the 
runway, as during taxii, take-off and landing, it is required to 
solve the vertical force equilibrium equations (2.39) and (2.40). 
Since, it is difficult to obtain a closed form solution for these 
coupled nonlinear differential equations in the general case, some 
numerical approach has to be adopted. For this a standard library 
routine using variable order, variable step size Adam’s method is 
selected. 

The initial conditions of state required for any numerical 
integration scheme to start, are known at the start of the run. 
At this instance the system is assumed to behave as a single 
degree of freedom system. When transition takes place to two 
degrees of freedom the number of variables increase and the 
appropriate initial conditions have to be set up to reinitiate the 
integration scheme. An iterative approach is developed for this 
purpose which is described later in this chapter. 

Besides this, a numerical scheme is also needed for the 
simulation of the mean and the random profile of the runway track 
as this information is required at each step for the integration 
of the system equations. 
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3.2 ln.it tal conditions for the tw degree case 

When "the aircraft comniences to move on the runway the 
modelled system is in the one degree of freedom case in which the 
shock struts do not telescope and both the sprung and the unsprung 
mass move as a single mass. The system functions as a single 
degree of freedom system until the forces in the shock strut as 
calculated by equations (2.47) and (2.48) exceed the preloads in 
the shock struts . 

The moment the preloads are exceeded the system starts 
behaving as a two degree of freedom system. At this point of 
time, it is necessary to establish the movements of the shock 
struts and consequently the displacements and velocities of the 
sprung and the unsprung mass. This can be posed as a problem of 
obtaining the root of a transcendental equation. An iterative 
scheme developed to obtain the displacements of the shock struts 
is presented below. 

Initially, a very small displacement is set in the strut, 
from which the strut telescoping velocity is obtained, based on 
the impact considerations discussed in the following section 3.3. 
Then these values of the displacement and velocity of the strut 
are used to find the force in the shock strut as shown in equation 
(2.35). The value of the force calculated is then compared to the 
force at which the system had transited from the one degree of 
freedom case to the two degree case. This process is repeated 
until the calculated force equals the force at transition. 

Various possibilities exists, sometimes the calculated force 
is less than the force at transition and then the values of the 
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stroke and velocity of the strut are increased and the reverse is 
true when the calculated force is larger than the force at 
transition. But it may also happen that the calculated force 
less but after some particular step it becomes larger 
than the force at transition. When this happens then linear 
interpolation is done to obtain the value of the stroke of the 
strut and the force is then calculated at that point. The process 
is repeated until the correct value of the stroke and velocity of 
the strut, corresponding to the force at transition is obtained. 

Once the strut telescoping velocities are determined then the 
angular velocity of the linkage is established and this is used to 
obtain the velocity of the differential motion of the two masses. 

Then using equation set (2.2) the initial conditions for the 
two degree case are established. 


3.3 Im.pact Consfi deration for determining the Strut Telescoping 
Velocity 

As seen in the previous section, a method has to be developed 
for obtaining the strut telescoping velocity because it is 
required to evaluate the force in the shock strut. 

Usually as the transition takes place very early during the 
motion of the aircraft and must happen in a very small interval of 
time, it could be considered to be occuring during an impact 
phase. During impact the shock struct force will vary with the 
strut movement causing a variable acceleration of strut closer. 

' Assuming the acceleration to change linearly with time 
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^ ^ (3.1) 

where is a constant to be determined and t is the time instance 
measured from the instance of the strut unlocking. 

Integrating equation (3.1) yields 


S = 




2 “ ^o 


where S = initial velocity 


- 0 (as it has just started from rest) 


Integrating once more 


(3.2) 


S = 


V 


c + S 
6 o 


(3.3) 


6 (S-S^) 


n 


(3.4) 


Substituting equation (3.4) in (3.2) gives 


3(S-S^) 
S = 


(3.5) 


Equation (3.5) is used to obtain the strut telescoping 
velocity during the iterative scheme, used for transition. 


3.4 Numerical Sinmlatiori. of thjs Random Profile 

The ground unevenness acts as a base excitation to the 
modelled system. The random part of the runway profile is 
simulated based on a method developed by Joel N Franklin C -(g ], 
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which confirms to the power spectral density of the runway track. 

The simulation is carried out in the following steps : 

First a sequence (x^) of independent samples from the uniform 
distribution on 0 < x < 1 is generated using a standard library 
routine. 

Using these, a sequence (co^) ©f independent sample from the 
Gaussian distribution with mean zero and unit variance is 
constructed by using the method of Box and Muller [20] . 

The sequence ig obtained from the formulas 


CO 

2n~l 


CO 

2n 


(-2. ^n Cos 2n x^^. n = 1, 2 

(-2. ^n ^2n’ ^= 1 , 2 .... 


^ (3.6) 


This one dimensional sequence . .. can be used to 
simulate samples from the n-dimensional Gaussian distribution by 
defining vectors 


CO) 

CO =: < 


CO 


CO 


n 


Cl) 


CO 


CO 


n-H 


CO 


2n 


(3.7) 


These will be used to determine the Gaussian random vector. 

Assuming that the given power spectral density $(n) can be 
represented with sufficient accuracy by a rational function, we 
can represent 5(0) as the ratio of two polynomials in Cl as 



2 
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i{Ci) = 


P(iQ) 

Q(iO) 


where the degree of P is one less than the degree of Q, 


The output x(t) is related to the input »(t) by 


= mW) • 


or equivalently, if D = d/dt 




The last equation is to be understood in the sense 
first solve the differential equation 

Q(D) . f (t) = <o(t) 

for the steady state solution $(t), and then obtain the 
signal 

x(t) = P(D).$(t) 


Let equation (3.11) in general have the form 


S''^(t) + a s'"' ‘'(t) + . . .+ a «{t) = co{t) (3 

1 n 

where <o(t) is white noise 

In order to compute equation (3.12) we will require 
of the state vector 


(3.8) 

(3.9) 

(3.10) 

that we 

(3.11) 
required 

(3.12) 

.13) 

samples 
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Z(t) = 


«(t) 

#’(t) 




(3.14) 


The vector 2 (t) satisfies the + v 

xies the stochastic differential 
equation ^-Lnerential 


where , 


<32(t) 

dt “ A.Z(t) + f(t) 


f(t) = 


0 

<*J(t) 


(3.15) 


(3.16) 


The vector 2(t) le a Gauselan Random Vector with a no .r, 
definite moment matrix =itlve 


M = E CZ(t). Z*(t)] 

The moment matrix M 

^ M can be computed without a 

the eigen values of A. The matrix m k 

axrix M has components 


(3.17) 


a knowledge of 
s of the fo 2 rm 


f “ 

l (-1 


> i+j odd 




«<j+i>/2}-i> * sven 


(3.18) 
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where the numbers m^, can be computed by solving the n 

linear equations 


(-1) - E 

k/2<q:<{n+k)/2 


(~1)^ a 


■n- 2 q+k ■ **q 


{ 0 , k=0, . . . ,n-2 

(3.19) 

, k=n-l 


where we define a = l 

o 


When M has been computed it can be factored by the Grout 
factorization 


M = TT 


(3.20) 


The required initial vector Z(0) - can be computed as 


Z(0) = 


(3.21) 


Suppose Z{t) has been computed for any t > 0, we will have to 
compute the succeeding state vector Z(t+At). 

X(t) = = exp At (3.22) 

be the solution of the initial -value problem 


dx(t) 

dt 

X(0) 


= A. X(t) 
= I 


(3.23) 


We need to know the values of the matrix exp (At) for the 
single value of time t = At. The solution of the differential 

equation gives the values of the matrix exp (At), at a particular 
time t = At 
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C = AM^ + A* 


(3.24) 


^ilex’s C is a positive semidef inite matrix 



0 

0 

6 


0 . . . . o' 
0 . . . . 0 

0 1 


(3.25) 


)nce [e T is known the Ljapunov equation (3.24) can be solved to 
obtain the moment matrix 

Having computed the matrix we factor it by the Grout 

ractor izat ion . 


We are now ready to compute a sample for r. 


Let t + At = v.At (3.26) 

) 

Let oi be one of the vectors defined in equation (3.7). We 
ihen compute 

Z(t+At) = Z(t) + r (3.27) 

If the polynomial P has the form 

P(p) = b .p“ + bp““* + ...+b^, b^ 5- 0, m < n (3.28) 

o 1 m o 

;hen the required scalar function x(t) with p.s.d, $(0) is 
!omputed as a linear combination 


x(t) = b Z (t) + b Z^(t) + ...+ b„ Z (t) (3.29) 

O TD lrl 1 Eft in 1 

Equation (3.29) is equivalent to equation (3.12) 


Chapter IV 


Results and Discussions 


4.1 1 rttradxLC t i on 

Response of a non-linear aircraft to the ground induced 
excitation has been considered during variable velocity runs 
landing and take-off. 

The various numerical data used for the system in generating 
the response values are provided in the Appendix. As a first 

study, the main shock strut is always considered locked and only 
the drag strut is allowed to telescope. 

The response obtained constitutes the mean values and the 
standard deviations for the displacement, velocity and 
acceleration of the sprung and the unsprung mass. Due to the 

physical constraints of time, the sample size used for obtaining 
the mean and standard deviation is restricted to fifteen in 
number. To smoothen out the irregularities in the response due to 
the small size of tlie ensemble, a ten point moving average is 


carried out over 

the generated 

response . 

Most 

of the times 

to 

obtain a more clear and vivid 

picture of 

the 

variations 

in 

the 

response for the 

presentation 

of graphs , 

the 

duration 

of 

the 


landing or take-off run is broken down into more than one interval 
of time. 

The duration of the landing run is taken as fifteen seconds 
and that of the take-off run as seventeen seconds. 
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For the landing case, variation in the values of the 
roughness constant of the runway track and the variations in the 
sink velocity of the aircraft at touchdown has been studied. For 
the take off runs the effect of variations in the roughness 
constant of the runway track or the response behaviour is studied. 

4.2 La.nd.tng Run with variation in Track. Roughness 

During the course of study the value of the sink velocity of 
the aircraft at touchdown was held constant at 350 cm/s. The 

value of the ground mean slope was also held constant at 1 in 
1000. The roughness costant C has been taken as 0.0001, 0.0005 
and 0.00002 in ft units. 

4.2.1 Mean Displacement of Unsprung mass 

The mean displacement response of the unsprung mass is 
presented in Figures 4.1a and 4,1b. The following observations 
are made: 

i) The track roughness has a very small effect on the mean 
response . 

ii) As time progresses during the landing run, the heaving 
of the unsprung mass ceases considerably beyond one 
second after touchdown and only a slight heave movement 
is discernible. 

iii) For the landing gear data used the aircraft does not 
'balloon’ after touchdown. 
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Fig. 4.1a Landing run. Roughness variation? Mean Displacement 

Response. 




Landing run, Roughnes-s variation; Mean Displacement 
Response. 
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4.2.2 Mean Displacement of the Sprung mass 

The mean displacement response of the sprung mass is shown in 
Figures 4.1c and 4. Id. Some observations are : 

i) The track roughness lias a very minor effect on the mean 
response. 

ii) At around 0.15 seconds after touchdown a resonance is 
observed. 

iii) Just after or very near to the time of impact, 
fluctuations are observed but with the passage of time 
these gradually decrease. 

iv) There is a steady increase in the response due to the 
motion of the aircraft over the mean slope. 

4.2.3 Standard Deviation of Displacement of Unsprung mass 


The 

standard deviation 

of 

displacement 

response of 

the 

unsprung 

mass is portrayed 

in 

Figures 4.2a 

and 4.2b. 

The 


♦ 

observations are : 


i) The influence of the variation in the track roughness is 
more easily noticeable than in the case of the mean 
displacement. 

ii) The graph for the time duration of 0 to 1 second shows 
two resonances occur ing at around 0.2 and 0.8 seconds 


from touchdown. 
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Fig. 4.1c Landing run. Roughness variation; Mean Displacement 

Response. 
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Fig. 4.1d Landing run. Roughness variation? Mean Displacement 

Response. 
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Fig. 4.Ea Landing run , Roughness v 

of Displacement Respons 
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Fig. 4.2b Landing run , Roughness variation? Standard Deviation 

of Displacement Response. 
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iii) The overall -trend of the graph shows a decreasing 
pattern of the standard deviation as time progresses 
. during the landing run. This holds true for the 
different roughness levels considered. At around 15 
seconds the value of standard deviation is as small as 1 
to 1 . 5 mm . 

4.2.4 Standard Deviation of Displacement of Sprung mass 

The response is shown in Figures 4.2c and 4. 2d. The 
inferences made are : 


i) 

As in the 

case of the unsprung mass. 

here 

also 

the 


effect of 

variations in the roughness 

levels 

of 

the 


track is ' 

easily seen. 




ii) 

Resonance 

is observed at around 0.35 

seconds 

and 

0.8 


seconds after touchdown. 

iii) At around 9.0 seconds also a large peah response is seen 
for a particular value of roughness constant, C = 
0.0001 ft. 

iv) The trend of standard deviation of the displacement 
shows a dropping pattern and at 15 seconds it has a 
value of about 1 to 1.5 mm. 

4.2.5 Mean Velocity of the Unsprung mass 

The mean velocity response of the unsprung mass is shown in 
Figures 4.3a and 4.3b. The salient points noted from the graphs 


are : 
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Landing run, Roughness variation? Standard Deviation 
of Displacement Response. 
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Fig. 4.3jb Lainding 'run. Roughness veiriationj Meein Velocity 

Response. 
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i) At different roughness levels the effect of the 

variations in the roughness levels of the track is not 
very significant, giving a maximvm value of the mean 
velocity of only about 2 cm/s. 

ii) After the initial impact phase, beyond 0.25 seconds the 
mean value of velocity is maintained at a very low 
value. The maximum value during the phase is 

approximately 2 cm/s. 

4.2.6 Mean Velocity of the Sprung mass 

The mean velocity of the sprung mass response is presented in 
Figures 4.3c and 4.. 3d. The noteworthy points observed are : 

i) As before in the case of unsprung mass, the effect of 

roughness variation is not appreciable. 

ii) After an initial drop in the mean velocity a second peak 
is observed at around 0.1 seconds from touchdown, 

4,2.7. Standard Deviation of the Velocity of the Unsprung mass 

The response is shown in Figures 4.4a and 4.4b. The points 
of interest observed are : 

i) Effect of roughness variation is more marked in the 
first one second of the run than during the later 


period of time. 
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Fig. 4.3c Landing run. Roughness variation; Mean Velocity 

Response. 
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Fig. 4.4a Landing run , Roughness variation?Standard Deviation 

of Velocity Response. 
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Landing run, Roughness var iat ion jStandard Deviation 
of Velocity Response. 
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ii) The broken up graph for one second duration shovfs the 

\ 

presence of five resonances at 0,025, 0.1, 0,2, 0.7 and 
0.85 seconds. 

iii) For a particular value of roughness constant, C = 0.0005 
ft., there is a resonance at around 8.5 seconds, 
otherwise there is no significantly large variation due 
to the roughness. The value of the standard deviation 
of velocity drops as time progresses. 

4.2.8. Standard Deviation of the Velocity of Sprung mass 

The standard deviation of the velocity response of sprung 
mass is presented in Figures 4.4c and 4.4d. The salient points 
are : 

i) Effect of roughness is more dominant for the first one 
second of the run. 

ii) For the first one second resonances are observed at 0.2, 
0 . 7 and 0 . 9 seconds . 

iii) For the C value of 0.0005 ft. a peak is observed at 8.5 
seconds . 

iv) The value also decreases continuously with time and at 
15 seconds is very close to 1 cm/s. 

4-2.9. Mean Acceleration of the Unsprung mass 

Figures 4.5a, 4.5b and 4.5c show the response of the mean 
acceleration of the unsprung mass. The main features of these 
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Fig. 4.4d Landing run, Roughness var iat ion jStandard Deviation 

of Velocity Response. 
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Fig. 4.5a Lan<iing run , Roughness variation jMean Acceleration 

Response. 
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Fig. 4.5b Landing run , Roughness variation ;Mean Acceleration 

Response. 
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figures are ; 

i) At the time of impact, the value of the mean 
acceleration is very high, roughly to the tune of 13 to 
14 g values upwards. The negative sign indicates that 
instantaneously the force due to this acceleration tries 
“to lift wheel off the track. 

ii) For the higher values of C, at C = 0.0005 ft and C = 
0.0001 ft, a strong decrease in acceleration takes place 
at around 0.03 second and then once again there is a more 
steady drop after 0.04 second. 

iii) Effect of the variation in roughness constant is not 
very dominantly observed. 

iv) The value of the mean acceleration keeps on reducing and 
at 15 seconds it is close to 0.01 g. 

4.2.10. Mean Acceleration of the Sprung mass 

Figures 4.5d to 4.5f show the response of mean acceleration 
of sprung mass. The interesting points observed are : 

i) Maximum value after impact is about 5.5 g, but it slowly 
keeps on reducing and subsequently steadies at 0.01 g. 

ii) There is some effect of the variation in the roughness 
level of the track. 

4.2.11. Standard Deviation of Acceleration of Onsprung mass 

The response of standard deviation of acceleration of the 
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Fig. 4.5d Landing run , Roughness variationjMean Acceleration 

Response. 





Fig, 4.5e Landing run , Roughness variationjMean Acceleration 

Response. 
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.5f Landing run, Roughness variationjMean Acceleration 

Response. 
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unsprung mass is shown in. Figuies 4.6a to 4.6c. The observations 
are : 

i) Upto a time of 0.04 seconds after touchdown the value is 
very high, , around 13 to 14 g. 

ii) The effect of roughness is also visibly observed upto 
0.04 second and from 0.6 to 1 second. 

iii) Resonance is observed at 0.6 and 0.95 seconds. 

iv) The value keeps on decreasing as time passes and from 1 

to 15 seconds there does not seem to be a strong 
influence of roughness level. 

4.2.12. Standard Deviation of Acceleration of Sprung mass 

Figures 4.6d to 4.€f show the standard deviation of 
acceleration response of the sprung mass. The points of interest 
are : 

i) Similar observations are made as for the unsprung mass 
with the only difference being the numerical values. 
The maximum value is about 1.2 g. 

ii) Resonance levels are also less than that for the 
unsprung mass. 

4.3 LandLing Rxirt with, ycariat ion irt Sink. Velocity 

For studying the effect of the variation of the sink velocity 
at touchdown, the value of the roughness constant was kept 
unchanged at 0.0001 ft. The mean slope is maintained unchanged at 
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Fig. 4.6a Landing run , Roughness var iat ion ^Standard Deviation 

of Acceleration Response. 
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Fig. 4.6b Landing runyRoughnessvariationjStandard Deviation 

of Acceleration Response. 
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Landing run , Roughness variation jStandard Deviation 
of Acceleration Response. 
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Fig. 4.6e Landing run, Roughness variation;Standard Deviation 

of Acceleration' Response . 
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Fig. 4.6f Landing run , Roughness variationjStandard Deviation 

of Acceleration Response. 
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1 in 1000. The values of the sink velooitv are taken as 350. 250 

and 100 cm/s. 

4.3.1 Mean Displaccmsn't of Unsprung mass 

The response of the mean displacement of the unsprung mass is 
portrayed in Figures 4.7a and 4.7b. The salient points worth 
noting are : 

i) The variation in the sink velocity has a marked effect 
on the maximum mean displacement of the unsprung mass. 
The effect is very clearly observed upto a time of 0.6 
second after touchdown. 

ii) As time progresses the heaving reduces considerably and 
after about 1 second the motion of the unsprung mass 
settles down and follows the mean slope, 

iii) Once touchdown is affected the aircraft does not get 
airborne and there is no balooning. 

4-3.2 Mean Displacement of the Sprung mass 

Figures 4.7c and 4.7d show the response of the mean 

displacement of the sprung mass. The points worth mentioning 
are : 

i) The effect of the variation in the sink velocity of the 
aircraft at touchdown is very clearly noticeable for a 
time duration of 0.7 second after touchdown. But later 
on this effect is not predominant. 
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Fig. 4.7a , Landing run, Sink Velocity variation? Mean 
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Fig. 4.7b Landing run, Sink Velocity variation; 

Displacement Response- 
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li) The mean displacement response reaches a maximum at 0.1 
sec after touchdown. Subsewentljr. it drops and shows a 
resonance at 0.2 sec. This behaviour is observed for 
all the three sink velocities considered. 

iii) At higher sink velocities of 350 and 250 cm/s the 
resonance value is smaller than that at the impact peak. 
The difference between these two grows smaller with 
reduction in the sink velocity. At the sink velocity of 

100 cm/s the resonance value becomes higher than the 
impact peak. 

Deviation of the Displacement of Unsprung mass 

The response of the standard deviation of the unsprung mass 
is shown in Figures 4.8a and 4.8b. The points observed are : 

i) The effect of the variation in the sink velocity of the 
aircraft is clearly visible. 

ii) Resonance of the displacement occurs very sharply at V 

s 

350 cm/s at around 0.8 second after touchdown. 

iii) The numerical value of the standard deviation of 
displacement shows a decreasing trend as time passes 
during the course of the run. This is true for all the 
sink velocities studied. 

Standard Deviation of the Displacement of Sprung Mass 
The response of the standard deviation of displacement of the 
sprung mass is shown in figures , 4.8c and 4.8d. Some salient 
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Fig. 4.8a Landing run, Sink Velocity variation? Standard 

Deviation of Displacement Response. 
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Landing run, Sink Velocity variation; Standard 

Deviation oP Displacement Response. 
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points are : 

i) The effect of the variation in the sink velocity is 
clearly brought out. As for the time duration of upto 
1 second after touchdown it is clearly seen that higher 
values of sink velocities give higher values of 
displacement, the same cannot be said confidently for 
the time duration of 1 to 15 seconds. 

ii) A sudden sharp peak value occurs for the sink velocity 
of 350 cm/s at 0.8 second after touchdown. Two more 
resonance peaks are observed at 9 seconds and at 12 
seconds . 

iii) The trend for the standard deviation is a decreasing one 
for all the sink velocities considered. 

4.3.5 Mean Velocity of Unsprung mass 

The response for the mean velocity of the unsprung mass is 
presented in Figures 4.9a to 4.9c. The points observed are : 

i) The variation in the sink velocity affects the mean 
velocity of the unsprung mass. 

ii) After a time duration of 0.3 second the mean velocity 
keeps on fluctuating about the zero level. 

4.3.6 Mean Velocity of the Sprung mass 

Figures 4.9d to 4.9f depict the response of the mean velocity 
of the sprung mass. The observations made are ; 
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Fig. 4.9a Landing run, Sink Velocity variation;Mean 

Velocity Response. 
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Landing run, Sink Velocity variationjMean 
Velocity Response. • 
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Fig. 4.9c Landing run, Sink Velocity var iationjMean 

Velocity Response. 
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Fig. 4.9d Landing run, Sink Velocity variation? Mean 

Velocity Response. 
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Fig. 4.9f Landing run, Sink Velocity var iat ion jMean 

Velocity Response. 
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i) The effect of variation in sink velocity is felt, but as 
in the case of the unsprung mass, a definite pattern in 
the response behaviour is not obvious. 

ii) After a time duration of 0.8 second after touchdown the 
mean velocity starts fluctuating at the aero level. 
This can be very well seen in Figure 3.9f. 


4-3.7 Standard Deviation of Velocity of Unsprung mass 

The response of the standard deviation of velocity of the 
unsprung mass is portrayed in Figures 4.10a to 4.10c. The 
essential points observed are : 

i) For the first 0.2 seconds of time duration the numerical 
values of the velocity of unsprung mass are higher for 
higher sink velocities and vice versa. 

ii) For the time duration of 0.2 to 1 second after 

touchdown for the sink velocity of 360 cm/s there are 
two reasonances at 0.7 and 0.9 seconds. For the sink 
velocity of 250 cm/s there is a single resonance at 0.95 
sec, whereas the sink velocity of 100 cm/s does not have 
a resonance. Thus the higher the sink velocities, the 
earlier is the possibility of the occurrence of 

resonance. 

iii) The trend of the numerical values is found to be a 

decreasing one. 
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Fig. 4.10b Landing run, Sink Velocity variationjStandard 

Deviation of Velocity Response. 
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Fig. 4.10e Landing run, Sink Velocity variationjStandard 

Deviation of Velocity Response. 



4.3.8 Standard Deviation of Velocity of Sprung mass 


Figures 4 . lOd and 4 . lOe present the standard deviation of the 
velocity response of the sprung mass. The observations made 
are • 

i) There is a clear effect of the variation of the sink 

velocity on the standard deviation of the velocity of 
the sprung mass. 

ii) For the time duration of 0 to 1 second resonance is 

observed at 0.2, 0.7 and 0.9 seconds for the sink 

velocity of 350 cm/s. For sink velocity V = 250 cm/s, 

s 

resonance is observed at 0.2 and 0.95 seconds, whereas 
there is no resonance for V = 100 cm/s. The 

observation is identical to that for the unsprung mass. 

4.3.9 Mean Acceleration of the Unsprung mass 

The response of the mean acceleration of the unsprung mass is 
shown in Figures 4.11a to 4.11c. The noteworthy points are : 

i) There is an effect of the variation in the sink 

velocity. Higher sink velocities give higher 
acceleration values at the time of impact. 

ii) The trend of the mean acceleration response is a sharply 
declining one. This is true for all the sink velocities 


considered. 
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Fig. 4.10d Landing run, Sink Velocity variationjStandard 

Deviation of Velocity Response. 
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Fig. 4.10e Landing run, Sink Velocity variationjStandard 

Deviation of Velocity Response. 




Fig. 4.11a Landing run, Sink Velocity' variationjMean 

Acceleration Response. 
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Fiq. 4.11b Landing run, Sink Velocity variationjMean 
' Acceleration Response. 
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Fig. 4.11c Landing run,-Sink Velocity variationjMean 

Acceleration Response. . 
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4.3.10 Mean Acceleration of the Sprung mas s 

The response of the mean acceleration of the sprung mass is 
presented in figures 4. lid to 4. Ilf. The observations made are : 

i) The effect of higher sink velocities, at a time just 
after touchdown is to increase the values of the 
accelerations of the sprung mass. 

ii) After 0.6 second the value of acceleration starts 
fluctuating around the zero level, 

4.3.11 Standard Deviation of Acceleration of Gnsprung naass 

Figures 4.12a to 4.12c present the response of the standard 
deviation of acceleration of the unsprung mass. The observations 
made are : 

i) The higher the value of the sink velocity, the higher 
the value of standard deviation of acceleration. 

ii) Resonance occurs for all the sink velocities at around 
0.6 and 0.96 seconds after touchdown. 

iii) The trend is a decreasing one. 

4*3.12 Standard Deviation of Acoeleration of the Sprung mass 

Figures 4.1 2d to 4.12f present the response of the standard 
deviation of acceleration of the sprung mass. The observations 
made are : 

i) The values of standard deviation of acceleration at 
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Fig. 4.11f Landing run, Sink Velocity variation;liean 

Acceleration Response. 















SD.DEV.OF AOCIELEPATION OF SPF=?UNO MASS 



12e Landing run, Sink Velocity variationjStandard 

Deviation of Acceleration Response. 
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Fig. 4.1Sf Landing run, Sink Velocity variationyStandard 

Deviation of Acceleration Response. 




higher sink velocities are higher than those at lower 
sink velocities. 

ii) The resonance characteristics are similar to those of 
the unsprung mass. 

iii) The overall numerical values are lower than those of the 
unsprung mass. 

1.4 Tahe-of/ run. with variation in Track koug-hness 

For the take-off run the variation in the roughness of the 


runway track is studied. 

The aircraft 

takes 

off 

from 

a track 

laving a constant mean 

slope. 

The 

values 

of 

the 

roughness 

jonstant C are taken as 

0.0001, 

0.0005 

and 0.00002 

in ft. 

, units. 


1.4.1 Mean Displacement of the Unsprung mass 

The response of the mean displacement of the unsprung mass is 
>resented in Figure 4,13a. The observations made are : 

i) There is very little or negligible effect of the 
variation in the roughness level of the track. 

ti) Initially at the start of the take-off run the unsprung 
mass keeps heaving but after about 2 seconds the heaving 
gives way to the gradual decrease in the tyre 

deflection. 

iii) The aircraft take-off takes place at around 15.5 seconds 
after the start of the run. 




Fiq. 4.13a Take-off run , Roughness variationjMean Displacement 


4.4.2 Mean Displacement of the Sprung mass 

Figure 4.13b shows the response of the mean displacement of 
the sprung mass. The observations are : 

i) Initially upto around 4 seconds there is heavy heaving 
of the sprung mass which dies down subsequently. Then 
there is a gradual increase in the response due to the 
motion of the aircraft on the mean slope. 

ii) There is a drop in the value of the response at around 

13 seconds indicating a valley after resonance. Then 
there is a sharp increase until the take-off point is 
reached . 

iii) The effect of the variation in the track roughness is 

negligibly small. 

4.4.3 Standard Deviation of Displacement of Unsprung mass 

The response of the standard deviation of displacement of the 
unsprung mass is represented by Figures 4.14a to 4.14c, The 
observations made from them are : 

i) Effect of variation of the track roughness is visible, 

but it is not possible to define a set pattern as a 
result of this variation. 

ii) The trend of the standard deviation response is an 

increasing one, but once the take-off is attained the 
standard deviation suddenly shoots up. 
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Fig. 4.13b Take-off run , Roughness variation?Mean Displacement 
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Fig. 4. 14c Take-off run, Roughness var iat ion ^Standard Deviation 

of Displacement Response - 
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Ul) A clear reaonance Is reallssd at around 0.6 seconds from 
the sstart of the run. 

standard Deviation of t he Dl3Place,.ent of 

Figures 4.14d to 4.14f show the standard deviation of 
displacement response of sprung mass. The following points can be 
highlighted : 

i) The effect of the variation of the track roughness is 
easily visible. Upto a time of 1 second, the increase 
in the track roughness leads to an increase in the 
numerical values of standard deviation. 

ii) A couple of resonances are observed, one at 0.6 second 
and the other at 1 second. 

A 

iii) The trend of standard deviation after the resonance is a 
decreasing one upto 12 seconds. After 12 seconds the 
deviation again starts increasing and once the take-off 
point is reached it immediately shoots up. 

Velocity of the Unsprung mass 

The response of the mean velocity of the unsprung mass is 
represented in Figures 4.15a to 4.15c. The prime observations 

made are : 

i) The effect of the variation in the roughness level of 
the trick is felt, but the roughness variation does not 
indicate a set pattern. 




Fig. 4.14d Take-off run , Roughness variation ;Standard Deviation 

of Displacement Response. 
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ii3 After the initial heavy fluctuations the mean velocity 
has very small fluctuations upto 12 seconds, but after 
that time the fluctuations again increase to a certain 
extent, 

iii) Once take-off point is attained the velocity drops off 
very sharply. 


4.4.6 Mean Velocity of the Sprung mass 

The response of the mean velocity of the sprung mass Is 
presented in Figure 4.15d. The observations made are : 


i) Effect of the variation of the track roughness is not 
very predominantly seen. 

ii) After about 2 seconds the mean velocity is seen 

fluctuating near the aero level. 

iii) Once the take-off point is reached there is a sharp 

decline in the velocity. 


4.4.7 Standard Deviation of Velocit y of Unsprung mass 

The response of the standard deviation of velocity of 
unsprung mass is shown in Figures 4.16a and 4.16b. The salient 
points noted are ; 

i) The effect of variation of track roughness is clearly 
visible in both the figures. 

ii) Two resonances are observed at around 0.5 and 0.75 
seconds after the start of the run. 



Fig. 4-15d Take— o^f run r Roughness variation? Meetn S/e loci ty 
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Fig. 4.16a Take-off run,RoughnessvarialiQn;Standard Deviation 

of Velocity Response. 





Fig. 4.16b Take-off run, Roughness variation yStandard Deviation 

of Velocity Response. 
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ill) After a time duration of 12 seconds there is one more 
resonance at around 14.5 seconds . 

4.4.8 Standard Deviation of Velocity of Sprung mass 

The response of the standard deviation of velocity of sprung 
mass is shown in Figures 4.16c to 4.16e. The interesting 
observations are : 

i) Effect of the variation of the track roughness is 
clearly visible, but it fails to indicate a set pattern 
of influence. 

ii) Resonances are observed at 0.5, 0.8 and 3.5 seconds for 
all the values of track roughness considered. 

iii) The trend of the standard deviation shows a decrease 

upto 12 seconds, but after 12 seconds there is again a 
large increase in the value. 

4.4.9 Mean Acceleration of the Unsprung mass 

The response of the mean acceleration of the unsprung mass is 
shown in Figures 4.17a and 4.17b. The noteworthy points are : 

i) The effect of the variation of track roughness is 

clearly visible after 0.3 second. 

ii) A resonance occurs at around 0.1 second after which 

there is a sharp decline in the mean acceleration. 
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Fig. 4.16ti Take-off run, Roughness variationyStandard Deviation 

of Velocity Response. 
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Fig. 4.17a Take-off run, Roughness variation ?Mean Acceleration 
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Fig. 4.17b Take-off run rRoughness variation;Mean Acceleration 
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iii) The fluctuations in the acceleration after 1 second are 
centered at the zero level, and at around 13 seconds one 
more resonance is observed. 

iv) After reaching the take-off point the value of the 
acceleration suddenly shows a sharp drop. 

4.4.10 Mean acceleration of Sprung mass 

Figures 4.17c and 4.17d show the response of mean 
acceleration of the sprung mass. The points observed are : 


i) 

The effect of 

variation 

of 

track roughness 

is 

not 

very 


predominant , 

only after 

12 

seconds some 

variation in 


response is 

observed 

at 

different values 

of 

track 


roughness. 

ii) The trend shows a fluctuation in the mean value upto a 
time duration of 3 seconds. Afterwards the value 
remains steady at the zero level upto 12 seconds. 

iii) Resonance is observed at around 13 seconds, but after 
that the value again starts oscillating about the zero 
level . 

iv) Once take-off is established the value shows a sharp 
decline. 

4 • 4 . 1 1 Standard Deviation of Acceleration of Unsprung mass 


Figures 4.18a amd 4-18b show the response of the standard 
deviation of acceleration of the unsprung mass. The observations 
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17c , Take~off run, Roughness variation;Mean Acceleration 

Response. 




MEAIJ ACCELEF?AriC'N OE K-tA- 



Fig. 4.17d Take-off run , Roughness variationjMean Acceleration 
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of Acceleration Response 




d.D.DE'/.OF ACCELER'i.TIQN OF UN5PRUM3 



Fig. 4.18b ‘ Take-off run , Roughness variationfStandard Deviation 

of Acceleration Response. 
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made are : 

i) Effect of variation, of the track roughness is very clear 
on the standard deviation of acceleration, 

ii) Resonances are observed at 0.45 and 0.55 seconds. 

iii) The trend of the graph shows a decreasing pattern after 
the later resonance has taken place. 

4.4.12 Standard Deviation of Acceleration of Sprung mass 

The response of the standard deviation of acceleration of 
sprung mass is shown in Figures 4.18c and 4.18d. The points 
observed are : 

i) The variation in the level of roughness of the track has 
some effect on the standard deviation of acceleration of 
the sprung mass. 

ii) During the first one second of the take-off run, three 
resonances at 0.2, 0.5 and 0,8 seconds are observed, 

iii) After 1 second the value of standard deviation for all 
roughness levels considered settles at around 0.01 g. 
This behaviour continues for about 12 seconds. Then at 
about 14 seconds after the start a resonance is observed 
and the value suddenly shoots up to 0.11 g. 

iv) After take-off the values are seen to drop sharply. 
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Fig. 4 . 18 c Take-off run , Roughness variation;Standard Deviation 

of Acceleration Response. 
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Fiq. 4.18d Take-off run, Roughness yar iationjStandard Deviation 

of Acceleration Response. 





















sink velocity 


u o 

(0 -H 
no 4 ^ 
CJ 

<d -H 
-p ^ 

CJQ 


(P tl0 

:j c 

iH O *H 
<0 «H no 
> 0 
^3 «5 
O iJ 
€ •Hi 
0 4^ tjo 
g <d 
•H ‘H 
X > 
to 0> 


•H 

u 

o 


S o TJ 


CS3 

o 

iH 

05 




















Table 4. 3 : Maximum values of Mean and Standard 


148 


o 

•H 

■P 

U 

a 

> 



u 

O 

0 

1 

0 

X 

<d 

o H 


-P bfi 
id d 
•H -H 

> h 

0 d 
p Td 


% 








o 

> 







M ^ 

w 

CS3 

to 

CO 

CD 

CO 

o 


p 

> 

• 

• 

> 

• 


< CO 

• 

00 

CO 

to 

CO 

CM 

tw 

« \ 

p 

t£> 

05 

to 

o 

CD 

o 

w 5: 

CO 

lO 

UO 


tH 


T-i 

p o 








W w 








o 

25 

CO 

CO 

CO 

to 

CM 


o 

< 

• 

• 

■ 

• 

• 

• 

< 

w 

tH 

CO 

o 

CO 

CM 

to 


s 

csa 

CO 

CM 

CO 

CO 

CO 



CO 

to 

CO 

xf 

1 

1 

1 


> 








m 


CO 

lO 

CO 

CO 

to 

£-* V— N 

p 

o 

CO 

o 


CO 

o 

M Cfl 

# 

CN3 

CO 

T-l 

CO 

CO 

r-4 

O ’V 

p 

• 


• 

• 

. 

• 

c s 

CO 

CO 

to 

CO 

lO 

to 

CO 

P u 








W 








> 


t-t 

CO 

CO 

CO 

CO 

CO 


< 

CO 

CO 


CD 

CO 




o 

CM 

C35 

CO 

CM 

CD 



CO 

1 

tw 

1 

CO 

1 

CO 

1 

C-* 

CO 

1 









2 









> 









CM 






w 

p 

lO 

C35 

o 

CD 

CD 

o 

o *--• 


CM 

CO 

o 

CO 

CO 

o 

< s 

p 



o 



o 

P o 

CO 

* 

• 

• 

• 

• 

• 

P 'W 


o 

to 


o 

to 

t- 

CO 








w 

22 







p 

< 

CM 

c- 

to 

CO 

CM 




CO 

to 

CO 

p 

CD 

CD 



lO 

to 

to 

■ 

• 

• 



. » 

• 

• 









CO 

CO 

CO 

CO 








CO 








w z 









* 






03 H e-* 



CM 

* 


CM 

dj CQ fa 

«rH 

to 

o 

r-l 

to 

O 

P 25 

O 

o 

o 

o 

O 

O 

O O 


O 

o 

o 

o 

o 

O 

cd o 


O 

o 

o 

o 

o 

O 



a 








2; 








D 



o 





« 



z 





P 

CO 


p 

CQ 




CO 

CO 


p 

CO 




% 

< 


p 

< 




D 

s 


CO 

S 




r-4 

d 

o 

IQ 

o 

0 

m 

to 


o 

-p 

d 

P 

* 

















Chapter V 


Conclusions 

In the present study the numerical data for the aircraft and 
the landing gear parameters contained in the appendix, are used to 
find the response of the aircraft. The aircraft is considered to 
be a rigid lumped mass with only heave degree of freedom and a 
single point input from the track. The response study is carried 
out for the variation in the roughness level of the track for the 
variable velocity runs - landing and take-off, and for the 
variation of the sink velocity during landing run. The response 
samples are then analysed to obtain the second order statistics of 
the response. The results and discussion for this study have been 
presented in Chapter - IV. 

5.1 Concluding Rornojrhs 

The salient conclusions drawn from the present study are 
listed below. 

i) Franklin’s method of numerical simulation of the random 
profile fits very well with the prescribed power 
spectral density (p.s.d.) of the track, for the entire 
frequency range. 

ii) The unsprung mass is subjected to much more severe 
forces than the sprung mass, as indicated by the higher 
values of accelerations for the unsprung mass. 

iii) During the study of the landing run, it is observed that 
the response values shoot up rapidly during the impact 
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phase, but with the passage of time, a decreasing trend 
is observed. 

iv) During landing, the effect of varying the roughness 

level of the runway track on the mean and standard 
deviation response is not very significant. 

v) Even for the small effect that roughness has on the 

response, a set pattern is not realised. The reason 
which can be attributed to this behaviour are the 

inherent nonlinearity of the system considered, 
nonstationarity of the track profile and the variable 
velocity runs of the aircraft. Only for the time 
duration of initial one second, in most cases, an 
increase in response with an increase in roughness level 
is observed. 

vi) During landing, the effect of the variation of sink 

velocity of the aircraft is predominantly observed in 
the heave displacements of the sprung and the unsprung 
masses. At higher values of sink velocities, the impact 
is more and as a consequence the masses are subjected to 
more heaving displacements. The acceleration response 
also indicates the same effect of sink velocity. Once 
the impact phase has passed, the definitive pattern of 
response due to sink velocity variation is not 
discernible. 

vii) During take-off the maximum values of mean and standard 

deviation respons^re lower than the values for landing. 


I y 1 


iii) During take-off, the effect of varying the roughness 
level of the track does not affect the mean response 
predominantly. The response standard deviation shows 
some effect, but a consistent pattern is not observable. 

ix) The response standard deviation of both the sprung and 
the unsprung masses show signs of resonance in many 
cases of landing and take-off runs. 


5.2 Suggestions for Futvire W'orh 

Extension of the present study could be carried out keeping 
in consideration one or more of the following, 

(1) Both the drag and and the main shock strut should be 
made to function simultaneously. 

(2) During the landing run, for the present study, lift was 
considered to be constant. Variation in the lift due to 
the variation in velocity must be considered. 

(3) Variations in the mean profile of the track can be 
studied , 

(4) Pilot's inputs can also be incorporated. 

(5) The ensemble size should be increased. 

(6) The differential equation solver should be such that 
during the process of discrete step-by-step solving, 
there is no back tracking, because this creates lot of 
computational complexities. 
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Besides these, extensions can be made in the actual system 
Bjodelling* Instead of the single point input from the track, 
multi point input can be studied. In that case the system model 
will ^ heave-pitch or a heave-pitch- roll model. Other 
variations may include unsymmetric landing, effect of side wind 
during landing and take-off, flexible aircraft structure (both the 
wing and the fuselage) and flexible runway pavement. 
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APPENDIX 

AIRCRAFT DATA 

I. GENERAL 


1. Gross weight of the aircraft 

2. Weight of the main wheels 

3. Wheel base 

4. Distance of aircraft c.g from 
main wheels 

5. Distance of aircraft c.g from 
nose wheel 

6 . Landing velocity of aircraft 

7. Take-off velocity of aircraft 


= 4154.0 Kg^ 

= 100.0 Kgj 

= 345.0 cm 

= 126.2 cm 

= 218.8 cm 
= 5400.0 cm/s 
= 5900.0 cm/s 


II.CHARECTERISTICS OF THE SHOCK STRDT OF THE MAIN GEAR 


1 . Coefficient of discharge 

2 . Coefficient of seal friction 

3. Mass density of fluid 
4. Orifice area 

5. Hydraulic area 

6. Pneumatic area 

7. Extended air pressure 

8. Extended air volume 

9. Index for air compression 


= 0.6 
= 0.01 

=0.85 gm/c.c 
= 2.405 sq.cm 
= 63.6 sq.cm 
= 49.5 sq.cm 
= 9.5 kg/cm^ 

= 1462.0 c.c 
= 1.3 



10. Initial oil volume 

11. Full stroke 

12 . Coefficient of friction at 
bearings 

13. Axial distance between bearings 

14. Coulomb friction 

15. Preload 


= 1950.0 c.c 
= 26.4 cm 

= 0.1 
= 20.0 cm 
= 0 ^ 

= 470.25 Kg^ 


III. MAIN TYRES 


1 . Radius 

2. Spin-up time 

3. Polar moment of inertia 


= 25.0 cm 
= 0.038 secs 
= 3217 kg-cm^ 


4. Load-deflection curve given 






